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Abstract 

We show that all bounded trajectories in the two dimensional classical 
system with the potential V(r, (p) = uj 2 r 2 + r - 2 ^ 2 ^ -f - are closed 

for all integer and rational values of k. The period is T — M- and does 
not depend on k. This agrees with our earlier conjecture suggesting that 
the quantum version of this system is superintegrable. 
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1 Introduction 



A recent article was devoted to an infinite family of quantum systems on a plane 
described by the Hamiltonian: 

H k = -d 2 - \ - \d% + + _ r ^— + -Jj- , (1.1) 
r r z r r z cos z kip r z sin kip 

where k is an arbitrary real number and (r,ip) are polar coordinates pQ. These 
systems were shown to be exactly solvable and integrable for all values of k. 
Integrability is assured by the existence of a second order integral of motion 



X ^ = Ll + ^— + ^—, L 3 = -id v (1.2) 
cos kip sin kp 



where L3 is the two-dimensional angular momentum. 

The existence of the integral Xj. implies the separation of variables in polar 
coordinates in the Schrodinger equation in quantum mechanics and also in the 
Hamilton-Jacobi equation in classical mechanics. 

Superintegrability in a two-dimensional system means that one further inte- 
gral of motion exists. Let us call it Yj.. 



It was conjectured in our previous article [T] that the system (1.1 1 is super- 
integrable for all integer values of k and that the additional integral of motion is 
of order 2k. This conjecture was proven for k — 1, 2, 3 and 4. If a = or /3 = 0, 
it was proven for k = 1 to 6 and 8. For k = 1,2 and 3, the system reduces 
to known superintegrable systems. For k = 1, this is one of the 4 systems found 
in [21 [3] , for k = 2 it was found in [3] and for k = 3 in [5] • 

In all cases (k = 1, 2, 3 and 4) the quantum Hamiltonian and the other quan- 
tum integrals of motion lie in the enveloping algebra of an underlying hidden 
Lie algebra [6] and this is presumably true for all integer values of k [1] . 

The purpose of this article is to analyze the corresponding classical system 
with Hamiltonian: 

? Pi 99 ak 2 Bk 2 ,„ „, 

r z cos z kp r z sin kip 

We show that all bounded trajectories for this system are closed and that the 
motion is periodic for all integer and rational values of k. 

We recall that in general, in an n-dimensional space, maximal superinte- 
grability means that the classical Hamiltonian allows 2n — 1 functionally inde- 
pendant integrals of motion (including the Hamiltonian) that are well defined 
functions on phase space. In classical mechanics maximal superintegrability 
implies that all bounded trajectories are closed and the motion is periodic [7]. 
The corresponding statement in quantum mechanics is that the existence of a 
non-abelian algebra of integrals of motion implies that the energy levels of the 
system are degenerate. In the two dimensional case that means that the energy 
depends on one quantum number, the wave function on two. This is indeed 
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the case for the Hamiltonian ( |1.1| for all integer values of k and is part of the 
justification of our conjecture PQ. 

In Section 2 we integrate the Hamilton-Jacobi equation for arbitrary values 
of k in polar coordinates, restricting to the case of bounded trajectories. In 
Section 3 we obtain an explicit expression for the bounded trajectories in terms 
of Chebyshev polynomials and show that they are all closed for integer and 
rational values of k. The period in time is T = For integer k with k > 2, up 
to k — 1 singular points can occur on the trajectories in with the tangent vector 
is not defined. The special case k = 1 is discussed in detail in Section 4. This 
is the only case in which the Hamilton-Jacobi equation separates in a second 
coordinate system, namely in cartesian coordinates. Section 5 is devoted to 
examples of trajectories illustrating all features presented in the earlier sections. 



2 Integration of the Hamilton-Jacobi equation 

2.1 The basic equations 

The Hamilton-Jacobi equation for the considered system in polar coordinates 
has the form 

fdS\ 2 1 fdS\ 2 2 2 1 / ak 2 (3k 2 \ „ „,„, 

\or J r z \oip J r z \cos z kip sm kip J 



We set the integral of motion (1.2) equal to a constant A. Since coordinates in 
(2.1 1 separate we look for a solution in which the classical action has the form 

S(r,<p,t)=S l (r)+S 2 (<p)-Et (2.2) 

and obtain 

r 2 f^) + u 2 r i -Er 2 = -A (2.3) 

Xt =(m' + ^ + j^ =A (2 . 4) 

\ dip J cos^ kip sin kip 

Since we are interested in the trajectories we do not need to solve for S\ and 
5*2, though that is obviously possible. Instead we follow the usual procedure [S] 
and calculate the derivatives of S with respect to E and A. We have 



V-A + Er 2 - u) 2 



i 



Si = / dr (2.5) 



, \J A cos 2 kip sin 2 kip — ak 2 sin 2 kip — f3k 2 cos 2 kip 

S 2 = : ; — : dip (2.6) 



and hence 





cos kip 


sin kip 


OS 


dSi 


t = 5! 






d~E 


~~b~E~ 




dS 




dS 2 


dA 


~~dA + 


dA ~ 



(2.7) 
(2.8) 
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where 5\ and 62 are constants. 



2.2 Integration of the radial part 

From d2.5h and \2.7\ we have 



t+Si = 



1 



r dr 



2J y/~A + Er 2 - u 2 r 4 



The roots of the denominator are 



(r U2 ) 2 = ^(E±y/E*-4w*A) 



Bounded motion T\ < r < r 2 is obtained for 



E 1 



4u 2 A > and A > 



Performing the integration in (2.9) we obtain 



—j (E + VE 2 - 4lj 2 A sin [4w(t + ^)] 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



Thus for w > 0, r 2 is indeed bounded and periodic with period T = £- for all 
values of a and j3 as long as (2.111 is satisfied. 

For ijj — the integral (2.9) must be evaluated differently and we obtain 

2 A + 4E 2 (t + S 1 ) 2 



E 



(2.13) 



We see that the trajectories for u> = are not bounded, independently of the 
values of alpha and beta. 

The limiting case A — corresponds to the harmonic oscillator for which all 
trajectories are bounded and periodic. In the other limiting case E 2 — 4lu 2 A = 

dS 

the integral (2.9) is not defined. Directly from (2.3) we obtain — — = 0, r 



dr 



1 E 



i.e. the trajectories are circles. 



2.3 Integration of the angular part 



From (2.51 and (2.6 1 we obtain 



dSi 
~dA 



dr 



rV-A + Er 2 - u 2 r 4 
sin kip cos ktp dip 



dS 2 = 

9 A 2 J yj A cos 2 kip sin 2 kip — ah 2 sin 2 kip — /3k 2 cos 2 kip 



(2.14) 
(2.15) 
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The integrals are readily evaluated to give 



dSi 
~dA 

dS 2 1 
dA 



I , -2A + Er 2 

-2A sin 2 kip + A-(a- f3)k 



aVa 



Ak^fA 



- 1 



^/(A-ak 2 + f3k 2 ) 2 -4k 2 f3A 



(2.16) 
(2.17) 

(2.18) 
(2.19) 

From the form of the Hamiltonian ( |1.3| we expect the trajectories to be re- 
stricted to sectors 

mr (n + \)n 



In addition to (2.111 we require 

(A - ak 2 + (3k 2 ) 2 - Ak 2 (3A > and (3 > 
A more symmetric way of writing this condition is 

[A-(a + f3)k 2 } 2 -4a(3k 2 >0 



< (p < 



n £ 



(2.20) 



2k r 2k 

The integral (2.15) was evaluated using the substitution z = sin 2 kip. The 
denominator in the integral vanishes for the values 



Zl,2 



(A - ak 2 + f3k 2 ) T \/P - ak 2 + (3k 2 ) 2 - AA(3k 2 
2A 



(2.21) 



Condition (2.18) is necessary to assure z\ < z < which is equivalent to 
(2.19). A final condition is < z\p. < 1 which implies 



A + k 2 (/3-a)>0 and a>0 



(2.22) 



3 The bounded trajectories 



3.1 Equation for trajectories for integer values of k 



So far we have established that bounded trajectories for the Hamiltonian (1.3 1 
are obtained if the constants of the system satisfy 



A>0, a>0, /3>0, uj > 0, A>k 2 \(3-a\, 
[A - (a + (3)k 2 ] 2 - Aa(3k 2 > 0, E 2 — 4uj 2 A > 



(3.1) 



Now let us obtain an equation for the trajectories. From (2.8), (2.16) and 
(2.17) we have 

k arcsin R + arcsin Uk = — / LkV~AS2 (3-2) 
where we have introduced 

R= - 2A + Er " and rr_-^sin 2 kp + A-(a-P)k 2 
r 2 y/E 2 - Auj 2 A y/[A- {a + (3)k 2 ] 2 - Aa(3k 2 
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for r 2 and z = sin 2 kip in the intervals 

-^(E - ^E 2 - Auj 2 A) < r 2 < -^(E + V ' E 2 - 4w 2 A) (3.4) 
2w z 

z\< z < z<i (3-5) 

where z 12 are as in (2.21). 
We have 

- 1 < R < 1, -1 < C/fe < 1 (3.6) 



Once R is known as a function of ip we obtain r(^>) as 
r((p) 



2A 



E - RVE 2 - 4lv 2 A 



We rewrite (3.2) as 

k arccos R = — arccos Uu + C (3.7) 

with 

C = ( k + 1 ^ n + 4kVAS 2 , < C < 2ir (3.8) 



Equation (3.7) can be converted into an algebraic equation for R in term of 
sin 2 k(p. We use the well-known formula [9] for the Chebyshev polynomials 



rp ( \ /, s (x + i-y/1 - x 2 ) k + {x- Wl - x 2 ) k 

T fe (x) = cos(fc arccos x) = (3.9) 

for k > 1 and obtain 
T k (R) = — — — = Uk cos C ± J 1 - f/£ sin C 



(3.10) 

We see that we can restrict the values of C to < C < it since the interval 
7r < C < 2ir will give the same trajectories. 

Let us write (3.10) as a polynomial equation with real coefficients for R (and 
thus for r 2 ). We do this separatly for even and odd values of k and obtain 



k = 2m 



T 2m = E ( 2 2 7) B- 1 )' Gj = ^ COS C ± V 1 - ^ C 



i=o v 7 j=0 
k — 2m + 1 : 

T2m+1 = 

'" ''2m ■ T ; 

, 2; 

Z=0 v 7 j=0 



E ( 2i ) E(-!) J ( • ) R2im - j)+1 = U *™ + l cos C±y/l- U 2 m+1 sin C 

(3.11) 
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The crucial point is that (3.11) does not contain any multivalued functions. 
It follows that R(ip) and ultimately r(p) is a periodic function with the same 

period r as sin 2 kip namely r — — and 

r^)=r(<p+^) (3.12) 
Hence all bounded trajectories are closed and the motion is periodic. 

7T 

We stress that r = — is the period of r as a function of ip. Both r and ip are 

7T 

periodic in time t with the same period, namely T = — (see (2.12)) and this 

period does not depend on k and coincides with that of the harmonic oscillator. 

Alternative (and equivalent) equations for the trajectories can be obtained. 
If we solve (3.10) for U k we obtain 

U k {p) = T k {R) cos C±y/1- T k (R) 2 sin C (3.13) 

Since we have \U k \ < 1 we can define r\ k by putting 

cosfe ± C) = U k (3.14) 

and rewrite (3.10) as 

T k (R) = cosr7 fc (3.15) 

In section 4 we shall plot the bounded trajectories for low values of k. The 
equation for the trajectories then simplify and we have 

k = l: R=UiCosC±^Jl - C/fsinC (3.16) 

k = 2: 2R 2 - 1 = U 2 cosC± ^1 - C/| sin C (3.17) 

fc = 3: 4i? 3 -3R=U 3 cos C ± ^1 - [7| sin C (3.18) 

fc = 4: 8i? 4 - 8i? 2 + 1 = C/ 4 cosC± ^1 - UfsinC (3.19) 

The trajectories are always closed curves within the sectors (2.20). They 
can never actually reach the boundaries of the sectors, unless we have a = or 
(3 = (or both) . The trajectories degenerate into line segments for U% = 1 (see 
(3.10)). The length of the segment is determined by the limits (3.4) and (3.5) 
which must be imposed. 

3.2 Singular points on the trajectories for k integer 

Let us consider a curve given in polar coordinates. Singular points on the curve 

dv 

are those where the derivative — is not defined. They satisfy the overdeter- 

dip 

mined system of equations: 

dF dF 
F{r,¥>) = 0, ^=0, — =0 (3.20) 



G 



In our case the equation of the curve is 



F k {r,<p) = T k (R) - U k cosCT\/l - U% = (3.21) 
and the derivatives satisfy 

dFom dTom dR 



dr dR dr 



=^( 2 2 7)d-^G)<"-^ 

1=0 v 7 j=o v/ 



dr 



dF2 m +\ _ dT 2m +i dR 
dr dR dr 



(3.22) 



f:( 2m 2 r)s<- i » , C) [2(m - i)+lli;2< ' 



1=0 v 7 3=0 

dF k _ cos Cy/l - Ul ± U k sin C dU k 



_ j} dR 
dr 



(3.23) 



From equation (3.3) we see that ^ vanishes only for r — > oo and only 

for </3 = — — , n = 0,±1,±2, .... These values correspond to points that do not 
2k 

lie on a trajectory. 

The singular points satisfy Uk = ± cos C that is 

. 2, A + fc 2 (/3 - q) ± + fc 2 (/3 - a)) 2 - AAfiW cos C 
sin K(^5 — — (3-24) 

and one of equations (3.23) for k even or odd respectively. 

QF k 

From (3.23) we see that — — is a polynomial of order k — 1 and will hence 

or 

have k — 1 zeros, all of them real. These zeros are actually singular points if 
the corresponding values of r and ip lie on the trajectories (see (3.21)). This in 
turn depends on the constants A, C and E i.e. on the initial conditions (for 
u>, a and (3 given). The statement thus is: the number no of singular points on 
the trajectories satisfies 

< n < k - 1 (3.25) 

For small values of k these formulas give: 
for k = 1 : 

dF-i dR , n 
-K- = 7* n <r <r 2 
or dr 

There are no singular points. 
For k = 2 : 

^=4^=0 
or dr 
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The possible singular point corresponds to r = 
For k = 3 : 

dF 3 , , dR 
The two possible singular points are R = ±^ i.e. 

r 2 - 



For k = A: 

^ = 16i?(2i? 2 - 1)^ = 
The three possible singular points are i?i = and i?2,3 = ±— p i-e. 



2A 2 2\/2A 



V2 



r 



3.3 Equation for trajectories for rational values of /c = m/n 

In the last subsection, we obtained the equation for the trajectories for integer 
values of k. We have established that those trajectories for the Hamiltonian 
(1.3) arc bounded if the constants of the system satisfy condtions (3.1). Now 
let us obtain the trajectories when k is rational. 

Here we set k — m/n for m,n integers. The bounded trajectories in this 
case are obtained under the same conditions on the constants of the system as 
those specified in the preceding subsection (3.1). 

Equation (3.7) becomes 

m arccos R— n arccos Uu + C (3.26) 
Using identity (3.9) we obtain 

T m (R) =T n {U k )cosC± ^/l-T n (U k ) 2 smC (3.27) 
Following the same arguments as in subsection 3.1, we see that R(<p) and 

7T T17T 

rUp) are periodic functions of the angle tp of period r = — = — . Again all 

k m 

bounded trajectories are closed and the motion is periodic. 

In Section 5 we will present the trajectories for the rational values k = 
1/2,1/3,3/2. 
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3.4 Comments on the range of the parameters 

The trajectories given by (3.10) or (3.11) (see also (3.16)-(3.19)) and (3.27) de- 
pend on three quantities E, A and C. E and A are values of the energy H^. 
(2.1) and the integral of motion Xk (2.4). The third integral Yk, responsible for 
superintegrability has so far not figured explicitly and is related to the constant 
C. The quantities A and E are restricted by relations (3.1) in order for the tra- 
jectories to be bounded. We allow C to take values < C < n. The trajectories 
could also be calculated differently using all three integrals of motion and just 
one of the Newton equations of motion. This may impose further conditions on 
the parameters involved. We discuss this in detail for fc = 1 in section 4. 



4 The special case k = 1 

The bounded trajectories for k = 1 are given by (3.16) with R(r) and Ui((p) 
defined in (3.3). The case k = 1 differs from that of all other values of k by the 
fact that the additional integral of motion Yf. for k = 1 is second order in the 
momenta and hence also leads to the separation of variables in the Hamilton- 
Jacobi equation. Explicitly we have |U [5] 



Y 1 = E x = Pi + lu 2 x 2 + -J (4.1) 



so Y± is the energy related to motion in the x direction (equally well we could 
have put Y\ = E<i = E — Ex or Y\ = E\ — E% where E% is the energy related to 
motion in y direction). The Hamilton- Jacobi equation in this case is multisep- 
arablc. Namely it separates in polar coordinates (for any k) and cartesian ones 
(for k = 1). For k = 1 it also separates in elliptic coordinates with an arbitrary 
focal distance of < / < oo. 

The equations of motion were solved in cartesian coordinates [3] and the 
solution corresponding to bounded trajectories in the notations of the present 
article are 



2uj 2 2cj 2 



x- = — + — \IE\ - 4aw 2 sin[4cj(t + C X )] (4.2) 



2 _ E 2 



V ' = 2uP + 2cW^ - W sm[M* + Ca)] (4.3) 
Above we have E\ + Ei = E and Ci, C2 are integration constants. From (|4.2|) 



and (4.3) we see that the trajectories are bounded and closed. They lie inside a 



rectangle 



Ei - s[E\ - Aau 2 ^ 2 ^ E x + s[E x - Aauo 2 



* X * ^2^ ^ 

E 2 - VE 2 - 4f3u 2 2 E 2 + ^E 2 - 40LJ* 

2lo 2 ~ V ~ 2lo 2 1 ' 

We see that 

2uj^fa < E x < E- 2cjV/3 (4-6) 
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Eliminating t from (4.2 1 and (4.3 ) we obtain an equation for the trajectories, 
namely 



V 



2 



E 2 ^E 2 - 4f3co 2 



2w 2 2uj 2 ^/E 2 -4aw 2 



cosS(~E 1 + 2lu 2 x 2 ) 

(4-7) 

± 2lj sin Sy/—a + E\x 2 — to 2 a 



for < S = Ci ~ C 2 < 7T. 

This equation is a fourth order polynomial equation involving the coordinates 
x and y. Viewing X clS el parameter for the curve we note that it must satisfy 
(4.5). The curve (4.7) has four branches, each lying entirely within one quadrant 
of the (cc, y) plane. For a > and (3 > it can never reach the coordinate axes 
(see (4.4) and (4.5)). 

In the limiting cases when sin S — the closed curve (4.7) collapses into a 
segment of a quadric. For S = this is a segment of a hyperbola and for S = n 
a segment of an ellipse (always inside the rectangle (4.4) and (4.5)). 

The trajectory (3.16) in polar coordinates can be made more explicit namely 

r 2 = ? 2A — — v (4.8) 

E - VE 2 -4lu 2 A(U cos C±Vl~U 2 sin CJ 

Again this is a fourth order polynomial curve (in r and sin ip) that degen- 
erates into a segment of a quadric for sin C = 0. Considering z — sin 2 ip as a 
parameter of the curve we must restrict to the interval (3.4). 

The trajectories given by (4.7) and (4.8) must coincide. This gives a relation- 
ship between cartesian constants {E\,E, S} and the polar constants {A, E, C}. 
We calculate r(t) and <p(t) from (4.2) and (4.3) and compare r(t) from (2.12) 
and tp(t) from (3.16). For given total energy E (and a,f3,uj) we obtain 



E X E 2 + 2{a + (3)lu 2 - y/(E 2 - Aau 2 ){E 2 - A[3uj 2 ) cos S 



and 

cos C = A(2E 1 - E) - E(a - j3) 

y/(A 2 + (a- (3) 2 - 2A(a + (3))(E 2 - Auj 2 A) 

The condition |cosC| < 1 is not necessarily automatically satisfied by the 
right hand side of (4.10) and imposes further conditions. First of all, the de- 
nominator in (4.10) must be the square root of a positive number which implies 

A>(V^+^P) 2 (4.11) 
Secondly, the equation cos 2 C = 1 is a quadratic polynomial in A and E\ 
uj 2 A 2 + A(E 2 -E 1 E~2(a+P)cj 2 ) + E 2 a+E 1 E(-a+(3) + (a-(3) 2 uj 2 =0 (4.12) 
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An example of the corresponding curve in the {E\ , A) plane is given on Fig. 
1. We have choosen E = 20, u> — 2, a — 3 and 0=2. The curve cos 2 C = 1 lies 
inside the rectangle 



E 



(Va+^) 2 <A< 4 ^ 2 

2u^fa <Ex<E- 2ay/uj 



We have < cos C < 1 in region I, — 1 < cos C < in II, cos C < — 1 III and 
IV, cos C > 1 in V and VI. 



2fl 



15 



10 




Rl 



10 



12 



14 



Figure 1: Region of bounded trajectories in (Ei, A) plane 



The dashed curve represents cosC = 0. At the intersection points of this 
curve and (4.12), we have cosC = g. 

5 Examples of trajectories 

In this section we present curves showing trajectories for specific choices of 
parameters u>, a and (3 and the integrals of motion E, A and C. We choose 
k= 1,2,3,4,1/2,1/3,3/2. 

From all of the figures we see that the trajectories do indeed lie within the 
sectors (2.20) and never reach the sector boundaries. From formulas in Sections 
2 and 3 and from the figures we see that we can actually restrict to one sector 
only e.g. the one with n = in (2.20). The trajectories in all others sectors are 
obtained by reflections in the sector boundaries. 
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The figures illustrate all features discussed in the previous sections. For 
k = 1, the trajectories have no singular points. The dependance on the constant 
C is quite spectacular. For all values of C the trajectories are fourth order 
curves. For C = 5 they resemble ellipses. For C — or C — ir the trajectories 
flatten out and finally degenerate into line segments (see Fig. 2-7). This happens 
for all values of k. For k = 2, (3.17) has two real roots. For each of them we 
present the trajectories in all 8 sectors for just one value of C (see Fig. 8 and 9). 
We see that the chosen trajectories have one singular point. For k = 3, (3.18) 
has three roots but only one is real. A typical set of trajectories is on Fig. 10. 
For k = 4, (3.19) has 4 real roots. We choose just one of them and present the 
trajectories in the first quadrant of the (x,y) plane (see Fig. 11). 

Trajectories for rational values of k, namely k = 1/2, 1/3,3/2 are given on 
Fig. 12, 13 and 14 respectively. They somewhat resemble Lissajous curves 
for an anisotropic oscillator. The symmetry with respect to reflection in the 
boundaries is obvious in Fig. 7-11. 

5.1 k=l 

We set E = 20, A = 12, u> = 2, a = 3 and (3 = 2. The trajectories are given on 
Fig. 2 to 7. 




Figure 4: C = f Figure 5: C = fj 
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Figure 7: Trajectories in the four 
sectors of the (x, y) plane for C = 

7T 

2 • 



5.2 k=2 



We set E = 50, A = 60, u) — 2, a — 3 and (3 — 2. The trajectories in the eight 
sectors for C = ^ are shown on Fig. 8 and 9. 



- 3 -2 

CP/ 
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1 2 3 
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-0.5 

/ 
/ - D5 




1.0 15 

a 




-1.0 

P 

-1.5 


; q 


\ 

\ 



Figure 8: Trajectories for the first Figure 9: Trajectories for the sec- 
root ond root 



5.3 k = 3 and 4 

For k = 3 we set E = 50, A = 100, w = 2, a = 3 and /3 = 2 and for fc = 4, 
£ = 16, A = 15, u = 2, a = | and /3 = i . The trajectories are shown on Fig. 
10 and 11. 
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\^ \ ' 








-3 -2 -1 


1 2 3 


a ' - 




/6., 





Figure 10: Trajectories for the 
real root for k = 3 




Figure 11: Trajectories for the 
first real root in the first quad- 
rant for k = 4 



5.4 k = 1/2 and 1/3 

In (3.28), we set E = 20, A = 24, uj = 2, a = 3 and f3 = 2. So for C* = 
trajectories are given on Fig. 12 and 13. There are no singular points. 
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5.5 *=? 

From (3.28), we set E = 15, A = 10, w = 2, a = \ and /3 = §. We obtain three 
roots for r 2 . Two of them are complex conjugate to each other. The trajectories 
are corresponding to the real root are for C = ^ are shown on Fig. 14. We 
observe two singular points. 




Figure 14: Trajectories for k = - 



6 Conclusion 

The present article, together with our previous one [1], provides a very strong 
indication that the Hamiltonian (1.3) is superintegrable for all integer values 
of k in classical and in quantum mechanics. In classical mechanics, the results 
of section 3.3 show that this is also true for rational values of k. Indeed we 
have shown in [T] that the system (1.1) is exactly solvable for all values of k 
thus supporting a previous conjecture |10j that all superintegrable systems are 
exactly solvable. A spectacular result is that the period T — ^- is not only 
independent of k but is the same as for the pure harmonic oscillator. 

We mention that the Hamiltonian (1.3) can be interpreted as representing a 
three body system on a line, once the motion of the center-of-mass is factored 
out. The same classical system was studied in |llj for oj — and a = 0. In 
this case the trajectories are never bounded (see (3.1)) and the motion cannot 
be periodic. A third integral of motion can still exist and the authors suggest a 
possible form of the additional integral [TT]. 

A rigorous proof of superintegrability requires the explicit construction of 
the third integral of motion for all integer and rational values of k and possibly 
even for any k > 0, k EM.. This is left for a future study. 
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